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Abstract 

We calculate the probability for rapidity gaps in the parton cascade 
for different approximations within the perturbative QCD and compare 
the results with recent measurements. The aim is to find out whether 
the dual connection between the parton and hadron final states - ob¬ 
served so far in various inclusive measurements - holds as well for the 
extreme kinematic configurations with colour sources separated by large 
rapidity gaps. A description of the data is possible indeed choosing the 
parameters of the cascade in the range suggested by recent analyses of 
the energy spectra (the k± cutoff Qo > QCD scale A ~ 250 MeV). 
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We calculate the probability for rapidity gaps in the parton cascade for different 
approximations within the perturbative QCD and compare the results with recent 
measurements. The aim is to find out whether the dual connection between the 
parton and hadron final states — observed so far in various inclusive measurements — 
holds as well for the extreme kinematic configurations with colour sources separated 
by large rapidity gaps. A description of the data is possible indeed choosing the 
parameters of the cascade in the range suggested by recent analyses of the energy 
spectra (the kj_ cutoff Qo ^ QCD scale A ~ 250 MeV). 


1 Introduction 


Recently the studies of pp collisions with high px jets at the Tevatron on 
and of dijet events in ep collsions at HERA 00 showed that a considerable 
fraction of these hard scattering events (of order 1-10%) contain large gaps 
in the rapidity distribution of particles. Such phenomena had been expected 
from the exchange of colour neutral objectsB. 

In e“'"e“ annihilation the primary parton process is the production of qq 
pairs followed by gluon bremsstrahlung and quark pair production. In this case 
the suppression of central gluon radiation would produce a rapidity gap for the 
partons but the recoiling qq like systems separated by the gap may become the 
source of nonperturbative hadron production to neutralize the color field; then 
no gap would remain between the observed hadrons. 

One way to generate large rapidity gaps between the final hadrons in e“'"e“ 
annihilation is through the production of colour singlet clusters separated by 
large rapidity gaps. Such configurations can be obtained perturbatively in hard 


+ e- 


qqqq 


processffiTeading to at least four partons in the final state, like e+e 
or qqgga\l\. In this case the hadrons are produced within each colour neutral 
qq or gg cluster and the large gap persists. As these processes involve a highly 
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virtual intermediate quark or gluon the rate for these large rapidity gap events 
is rather small though. One finds that the rate keeps decreasing with increasing 
gap size, contrary to the case of pp or ep collision. 

Such an unlimited decrease of the gap probability has been seen already 
some time ago§. The recent analysis at the resonance by the SLD Collabo¬ 
ration!^ has shown the fall of the gap probability over five orders of magnitude 
(see also the first results from ALEPHt^i). At the same time the observed 
absolute size of the gap probability exceeds the expectations from the above 
calculations Bill by about two orders of magnitude. This observation excludes 
this type of hard colour neutralization to be the dominant mechanism for the 
production of large rapidity gaps. 

In this paper we investigate how the fraction of rapidity gap events f{Ay) 
for the parton cascade (without hadronization) compares with the experimen¬ 
tal results on hadrons. At first sight one might expect that the decrease of 
the gap fraction for hadrons and therefore in the experimental data is steeper 
than the one of the parton cascade because of the additional suppression from 
removing the hadrons emitted into the gap during the hadronization phase. 
We will see however that this effect strongly depends on the scale Qo for which 
the parton cascade is terminated. 

There are many phenomena in multiparticle production which are in favour 
of a ver y s imple hadronization picture, called local parton hadron duality 
(LPHD) 1111 , which compares the observed sufficiently inclusive particle den¬ 
sities and correlations directly with the calculations at the parton level, as¬ 
suming in effect the pr oduction of a hadron close to the parton in phase space 
(for reviews, seer^tru). In this approach the parton cascade is evolved down 
to the low scales Qo ^ ^ of a few hundred MeV; Qo denotes here the cutoff 
in the transverse momentum of the emitted partons and A the QCD scale for 
the running coupling. 

This calculation of the gap rate is similar in spirit to the “clan” model 
analysis within a simplified parton shower interpretationllll the primary inde¬ 
pendent emission of hadronic objects leads to a Poisson distribution with the 
gap rate related to the probability for no emission by an exponential. 

In the following we calculate the gap probability analytically in different 
approximation schemes of perturbative QCD and compare to the data for the 
previously determined values of Qo and A. This will teach us whether the 
parton hadron duality works in the more extreme, quasi exclusive configura¬ 
tions with large rapidity gaps as well; in this case a confirmation of the duality 
approach is even less expected than in case of inclusive observables. 
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2 Rapidity gaps in the parton cascade 

A gap in the rapidity distribution occurs if the initial qq pair does not ra¬ 
diate into the respective rapidity interval. The probability for no radiation 
into a certain angular interval is given in field theory by the exponential Su- 
dakov form-factor lU, originally derived in QED. A problem very similar to 
the rapidity gaps considered here occurs in the calculat ion erf multijet rates in 
e+e“ annihilation applying the Durham/fcj^ algorithmEjO the 2-jet rate is 
again given by the Sudakov form factor which represents the probability for 
no-parton-production above a certain cutoff (resolution parameter). 

In our application we consider the rapidity gap having no gluons inside 
the angular or rapidity interval above the transverse momentum cutoff Qo 
which corresponds to a hadronic scale of a few hundred MeV as outlined in 
the introduction. 

Let us consider specifically the angular interval between 0i and 02 (0i < 
02) where 0 is measured with respect to the quark direction. The rapidity 
is then obtained from y = — Intg^. Let us further denote the probability for 
emission of a gluon at an angle 0' with the energy lo' off a parent parton A 
(either a gluon(g) or a quark(q)) as pa(w',0') = dnA/dio'dQ'. 

Then, the Sudakov form factor for the angular ordered cascade is given by 

A_4(T’, 0,Qo) = exp(-WA(T’, 0 ,Qo)) (1) 

wa{P,Q,Qq) = [ dio' f dQ'pAiuj',0'), ( 2 ) 

J Jkj_ OQo 

and this represents the probability that no gluon is emitted within the cone of 
half angle 0 from the parent parton with the energy P = Qj2 with transverse 
momentum above Qo- In particular, 

Aa(02)/Aa(0i) = exp (-^^(02) + wa(0i)) (3) 

represents the probabilty that there is no emission of a gluon with emission 
angle between 0i and 02 . 

These rates have been calculated in different approximations which we dis¬ 
cuss in the following. 

The double logarithmic approximation (DLA) 

The simplest approximation takes into account only the leading contributions 
from the angle and energy singularities of the gluon emission. In this case the 
total number of gluons radiated from a primary parton into a forward cone of 
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half angle 0 as in (|^ is given in a small angle approximation by 

Ca dLo' dQ’ ^ 

wAMM = — eT^of^ e) 

or, using logarithmic variables 

Ca n 

wa{Y,X) = — drj dr]'-/Uv') 


with 


r = ln5^, A = ln.^° 


Qo ’ A 

and the DLA anomalous dimension 7 g = 2Ncas/'!T for running coupling a. 




4iVc 


'foiv) = ittT’ = -T^’ ^ - ^nf 


11 

y 


(4) 

(5) 

( 6 ) 
S 

(7) 


r] + X b 

where Nc and n/ denote the numbers of colours and flavours respectively. 
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Furthermore Ca = | for quarks and Ca = Nq = 3 for gluons. One finds 


wa(Y,X) = + X)ln^^ 


-Y} 


( 8 ) 


The probability for a gap without primary gluons is then obtained as 

/a(0i, 02) = (9) 


for two arbitrary angles 0i <02 < § measured with respect to the quark 
direction. In this approximation the jets evolve independently in both hemi¬ 
spheres and the respective probabilities factorize. In particular, in the case of 
a symmetric angular interval (0g, — 0g) the gap probability is given by 

fAiOc) = ( 10 ) 


In the present small angle approximation the rapidity is related to the 
angle by y !« — In ^; in case of the symmetric rapidity interval of full width 
Ay, the relevant gap angle is 0 g = 2exp(—Ay/2) which yields 0 g = 2 for 
the limiting case Ay = 0. The use of the small angle approximation for such 
large angles is formally not allowed, but we will see below that the effect is 
numerically small (radiation occurs dominantly at small angles). The gap 
probability /a (Ay) following from ( 0 ) corresponds to an almost exponential 
decrease for the relevant parameter range. 
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It is interesting to note that for fixed coupling there would be a flattening 
of the gap fraction /^(0 g) for large gaps Ay as can be derived from 

wa(Y) = (fixedos) (11) 

This difference comes from the strong radiation near the kinematical limit 
k± ~ Qo ^ A in case of running as which becomes important for small lim¬ 
iting angles 0. In other words, for large gaps the strong radiation near the 
kinematic limit is kinematically allowed and has to be suppressed which makes 
the gap probability much smaller than in the case of fixed as where this radi¬ 
ation is not equally strong. 


The modified leading logarithmic approximation (MLLA) 

Besides the double logarithmic terms one also takes into account in this approx¬ 
imation the leading single logarithmic corrections. In this subsection we first 
derive the full 0(as) emission probability for e“'"e“ annihilation from which 
the various approximations beyond DLA can be derived. 

We denote the energy fractions of q,q and g as x+,x- and Xg and take 
the jet (thrust) axis along the (say, left-moving) q. The energy distribution is 
then given by 

(see, for example. Ref. Ill where also useful kinematic relations are given). 

In order to evaluate the gluon emission probablity from a parent quark 


Pg(w',0') = 


d{x+,x-) 

d{uj',Q') 


dn 

dx+dx- 


(13) 


we need to work out the Jacobian factor to move from {x^,x-) to (uj',0') 
using the kinematical relations 


x+ 


1 — Xg(l — Xg/2)(1 + COS0') 
1 — Xg(l + COS0')/2 
_ 1 ~ _ 

1 — Xg(l + cos 0')/2 


(14) 


with Xg = 2uj'IQ. This can be done ecomically by noting a;_ = 2 — a:+ — Xg 
and we obtain a simple formula 


J = 


d{x+,x_) 

5(w',0') 


2 f dx-\ a;^XgSin0' 


(15) 
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The final form of pq{i-o', 0') from parent quark is 


/ / C'f ,a;?_x„sin0' + xt 

® Tvl-TTI- ul-V' 

4-/Vc’ Q{l-Xg) (1 - a:+)(l - x-j 


(16) 


Using the kinematical relation x^ = 4(1 — a;+)(l — x_)(l — Xg)jx'^ one finds 
the emission density 


Pg(w',0') = — 


1 


uj' sin 0' Nc 


Cf 2 

7o -^- 


(17) 


From the exact result we can easily obtain the approximation for small 
gluon transverse energy fraction xt, i-C. for a;_ —> 1, 


Cf 


1 


2a;' ',2 


Pq(w',0') = -^70 , ■ n' 
Nc ^ Sin 0' 


i + (i-^) 


(18) 


1 _!_ 2 

This formula contains the DGLAP splitting function as 

factor and can be seen to correspond to the expression 


dn 


as 1 


Pqq{x+), 


dx+dXT TT XT 
for, with xt = 2u)' sinQ'/Q^ the appropriate Jacobian factor is 


(19) 


i9(a;+,a;T) _ 2(1 - x+) 


( 20 ) 


The Eq. ( |l8|) differs from the DLA result by the inclusion of the non¬ 
singular DGLAP term and the correct angular dependence. We have computed 
Wq by integrating the density as in (||) numerically with the bounds 


Qo 

sin 0 


<uj' < P, 


sin-1 ^ < 0/ < 0 
UJ' 




( 21 ) 


This yields an improved result for the gap probability calculated again as in 
(i) or (0) which we refer to as “MLLA-n”. 

An analytical approximation (“MLLA-a”) can be obtained by simplifying 
the integral over the nonsingular parts of the splitting function^. One obtains 


..DLA 


In - e^{Ai(A) - Ai(F + A)} 


( 22 ) 
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where Ei{z) = dte ^ jt and denotes the expression in Eq. (||). The 

results from these calculations will be discussed in the next subsection. 


Comparison of analytical results and a Monte Carlo calculation 
The result of our computations in the different approximations are displayed in 
Fig. 1. For the parameter A = 0.1 at A = 0.244 GeV we show the DLA result for 
the symmetric gap and the MLLA results in the analytic approximation 
and in the numerical evaluation of (||) with (18). Also shown as data 
points are the results from the ARIADNE Monte Carlo at the parton level 111 
for the same parameters. This Monte-Carlo generates the quark gluon cascade 
as a sequence of dipole emissions and it is terminated as in the analytical 
calculations by a transverse momentum cutoff El. Furthermore, there is no 
limitation to the relative magnitude of the parameters (besides Qo > A). 

As can be seen from the figure the MLLA corrections are not very large: 
about a factor or two after a decrease of the gap fraction by 5 orders of magni¬ 
tude. Of the same order is the deviation from the parton Monte Carlo which 
takes into account effects not considered here: the correct angular recoil from 
the primary gluon emission also for large angles and the spillover of secondary 
gluons into the gap if the primary gluon is produced ouside the gap. This 
effect could be responsible for the difference between the MC and the more ac¬ 
curate MLLA calculations. Surprisingly, the simple DLA and the MC results 
are quite close to each other. 

The figure also shows the very strong dependence on the cutoff Qo for 
fixed A. Indeed, decreasing the transverse momentum cutoff from Qo = 1 GeV 
down to Qo = 0.27 the gap rate at Ay = 6 decreases by 4 orders of magnitude. 
This comes from the singular emission probability dk±/kj_as{kj_) enhanced by 
the running a^: the closer we come to the pole at fcj_ = A the more difficult it 
is to avoid the gluon bremsstrahlung. 


3 Comparison with data 

Apparently the prediction for the gap probability in the parton cascade de¬ 
pends sensitively on the parameter Qq. A determination of both parameters 
A and Qo within the MLI^-LPHD approach has been perfomied in the study 
of particle energy spectra cll applying the moment analysis E3. This study has 
shown that the two parameters are very close to each other and there is an 
upper limit to the A-parameter: 


A < 0.1, 


(23) 


7 




Figure 1: Fraction of rapidity gap events in e'*"e“ annihila tion as a function of the full 
width of the symmetric gap. The DLA predictions from eq. ( [l0[ ) for different values of the 
transverse momentum cutoff Qo? s,nd for one value of Qq also the analytical (a) and the 
numerical (n) MLLA results (see text). The data points refer to the calculations of the 
ARIADNE Monte Carlo at the parton level for the corresponding parameters. 




whereas the absolute scale was found to be Qo — 270 MeV. A lower limit on 
A could in principle be obtained from the study of the mean multiplicity in 
a calculation of yet higher accuracy. As central value for our prediction we 
choose therefore Qo= 270 MeV and A= 0.1 (or A = 0.244 GeV) but smaller 
values of A are allowed as well. 

In Fig. 2 we compare such results from the DLA (which agree closely with 
those from the parton MC) with the experimental measurement of the fraction 
of rapidity gap events from the SLD Collaboration 11. These data also contain 
a contribution from the r-lepton events which is important for the large gaps 
Ay. The data are well described by the JETSET Monte Carlo El including 
the r-decays. The Monte Carlo predictions without r-decays are also shown 
in the figure and indicate a steeper falloff for large Ay. Furthermore it should 
be noted that the data refer to rapidity gaps between charged particles. From 
the Monte Carlo one would expect a sightly steeper dependence if the gap is 
calculated from all final state particlesc!l 

One can conclude from the figure that the parton cascade with the low 
values for the scale parameters Qq and A obtained from previous analyses 
gives already enough suppression for large rapidity gaps without any additional 
hadronization; values A = 0.05-0.1 provide an adequate description of the gap 
distribution (our predictions refer to gaps in the final state of all charged and 
neutral particles). 


4 Discussion and further predictions 

The data are consistent with a description of the hadronic final state in terms 
of a parton cascade terminated at a rather low cutoff of the order of the final 
state hadron masses without any additional hadronization phase. This could 
be interpreted in favour of a dual correspondence of parton and hadron final 
states and a soft colour neutralization mechanism . 

A go od description of the data has been obtainedH by the JETSET Monte 
Carlo In this model the parton cascade is terminated at a cutoff Qo 
0(1 GeV) and followed by a hadronization phase. According to Fig. 1 the 
suppression of large rapidity gaps at a scale Qo ^ 1 GeV is not so strong and 
most of the suppression must come from the resonance and particle production 
in the hadronization phase. In our parton calculation the strong suppression 
of large gaps again comes mainly from the last stage of the cascade evolution 
where the radiation becomes stronger because of the infrared enhancement and 
the increasing coupling constant as. It appears that the average properties of 
the hadronization phase can be well represented by the parton bremsstrahlung 
cascade with running coupling. 
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Figure 2: Fraction of rapidity gap events in annihilation as a function of the full width 

of the symmetric gap. The data points refer to the measurement of gaps between charged 
particles (r-lepton events included) by jthe SLD Collaboration □. The dashed histogram 
shows the expectation from the JETSETc3 Monte Carlo without r-leptons. Also shown are 
the DLA predictions for the gaps of a parton cascade as in Fig. 1. 
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In the following we note a few further predictions from the perturbative 
approach to the gap rate. 

1. There is a difference in the gap fraction of quark and gluon jets according 
to the respective colour factors in (||). This can be tested, for example, 
by studying gaps in jets at high pT in pp collisions between angles of, 
say, 01 < 02 < 45° using eq. (^). The slope roughly doubles when going 
from quark to gluon jets. This measurement would test the connection 
of the gap probability to single gluon emission. 

2. In the approximation considered the evolution of the parton cascades in 
both hemispheres is independent, so the corresponding rates factorize. 
For example the slope should double when going from the gap (0 g, f) 
to (0G,7r - 0g). 

3. The eq. (^) also predicts the energy dependence of the slope. With rising 
jet energy the gap distribution gets slightly steeper: at (5=200 GeV the 
slope is larger by 4 % in comparison to (5=90 GeV. 

It will be interesting to study the gap events further in order to test the pro¬ 
posed perturbative approach with soft colour neutralization. 
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